Correlation functions and excitation spectrum of the frustrated ferromagnetic spin-| 

chain in an external magnetic field 
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Magnetic field effects on the one- dimensional frustrated ferromagnetic chain are studied by means 
of effective field theory approaches in combination with numerical calculations utilizing Lanczos 
diagonalization and the density matrix renormalization group method. The nature of the ground 
state is shown to change from a spin-density-wave region to a nematic-like one upon approaching 
the saturation magnetization. The excitation spectrum is analyzed and the behavior of the single 
spin-flip excitation gap is studied in detail, including the emergent finite-size corrections. 



I. INTRODUCTION 

The interest in helical and chiral phases of frustrated 
low-dimensional quantum magnets has been triggered by 
recent experimental results. While many copper-oxide 
based materials predominantly realize antiferromagnetic 
exchange interactions, several candidate materials with 
magnetic properties believed to be described by frus- 
trated ferromagnetic chains have been identifiedji'^' 3 -'^'^ 
including Rb 2 Cu 2 Mo 3 Oi2 (Ref. H), LiCuV0 4 (Refs. i, 
HUH), and Li 2 ZrCu0 4 (Ref. i). The frustrated anti- 
ferromagnetic chain is well-studied)£ but the magnetic 
phase diagram of the model with ferromagnetic nearest- 
neighbor interactions remains a subject of active theoret- 
ical investigations. 8 ' 9 ' 10 ! 11 

In this work we consider a parameter regime that is 
in particular relevant for the low-energy properties of 
LiCuV04, corresponding to a ratio of J\ ~ — 0.3 J 2 be- 
tween the nearest neighbor interaction J\ and the frus- 
trating next-nearest neighbor interaction J 2 > 0. As 
the interchain couplings for this material are an order 
of magnitude smaller than the intrachain ones; 3 - we an- 
alyze a purely one-dimensional (ID) model. Apart from 
mean-field based predictions^ the nature of the ground 
state in a magnetic field h is not yet completely known. 
Therefore, combining the bosonization technique with a 
numerical analysis we determine ground-state properties 
and discuss the model's elementary excitations. 

The Hamiltonian for our ID model reads: 



L 

H = J2(JiS x - S x+1 + J 2 S X - S x+2 ) -hJ2s z x , (1) 

x=l 



where S x represents a spin one-half operator at site x. 

Bosonization has turned out to be the appropriate lan- 
guage for describing the regime |Ji| <C J 2 of Eq. fT}. 
This result has been established by studying the magne- 
tization process yielding a good agreement between field 
theory and numerical data* 9 - The derivation of the ef- 
fective field theory is summarized in Sec. HH Here, we 



extend on such comparison of analytical and numerical 
results and further confirm the predictions of field the- 
ory by analyzing several correlation functions in Sec. IIHI 
Then, in Sec. IIV1 we numerically compute the one- and 
two-spin flip excitation gaps and compare them to field- 
theory predictions. Finally, Sec. IVl contains a summary 
and a discussion of our results. 



II. EFFECTIVE FIELD THEORY 

We start from an effective field theory describing the 
long- wavelength fluctuations of Eq. (fT]). In the limit of 
strong next-nearest neighbor interactions J 2 3> | Ji|, the 
spin operators can be expressed as: 



S z a (r) 
S~(r) 



m + c(m) sin \2kpr + V^iripa} +••• 

(-iy G - ie "^ + ■■■ . (2) 



kp = (i — m)7r is the Fermi- wave vector and a — 1, 2 enu- 
merates the two chains of the zig-zag ladder. In relation 
with Eq. (H|), note that Si(r) = S( x+ i)/ 2 (S 2 (r) = S x / 2 ) 
for x odd (even). <f> a and a are compactified quantum 
fields describing the out-of-plane and in-plane angles of 
fluctuating spins obeying Gaussian Hamiltonians: 



n= v - J dx{^(d x ^ a ) 2 + K(d x e a ) 2 }, 



(3) 



with \4> a (x), 9 a {y)] = iQ{y — x), where O(x) is the 
Heaviside function. Sub-leading terms are suppressed in 
Eq. m is the magnetization of decoupled chains, re- 
lated to the real magnetization M of the zig-zag system 
by: 



M ~ m 1 



2K(m)Ji 
TTv(m) 



(4) 



K(m) and v(m) are the Luttingcr liquid (LL) parame- 
ter and the spin- wave velocity of the decoupled chains, 



2 



respectively. The nonuniversal amplitude c(m) appear- 
ing in the bosonization formulas ^ has been determined 
from density matrix renormalization group (DMRG) 
calculations^ Note that in our notation M — 1/2 at 
saturation. 

Now we perturbatively add the interchain coupling 
term to two decoupled chains, each of which is described 
by an effective Hamiltonian of the form Eq. ([3]) and fields 
4>i and 6i, i = 1,2. For convenience, we transform to 
the symmetric and antisymmetric combinations of the 
bosonic fields tf>± = ((f>i±fa)/V2 and 9± = (9 1 ±6 2 )/V2. 
In this basis and apart from terms of the form ([3]), 
the effective Hamiltonian describing low-energy proper- 
ties of Eq. (Q} contains a single relevant interaction term 
with the bare coupling g\ oc J\ <C v: 



H e s = Hq +H + gi J dx cos (hp 



(5) 



and the rcnormalizcd LL parameters K± are, in the weak 
coupling limit: 



K ; = K ( IT Ji — 

IT V 



(6) 



K + is the Luttinger-liquid parameter of the soft mode 
of the zig-zag ladder. The Hamiltonian ((SJ) represents 
the minimal effective low-energy field theory describing 
the region J2 ^> | J\\ of the frustrated FM spin- 1/2 chain 
for M ^ The relevant interaction term cos y/8n<j>- 

opens a gap in the </>_ sector. Since S* +1 — S* ~ d x <fi-, 
relative fluctuations of the two chains are locked. This 
implies that single-spin flips are gapped with a sine- 
Gordon gap in the sector describing relative spin fluctua- 
tions of the two-chain system.— Gapless excitations come 
from the AS Z = 2 channel, i.e. only those excitations 
are soft where spins simultaneously flip on both chains. 
DMRG results show that this picture applies to a large 
part of the magnetic phase diagram;* 



III. CORRELATION FUNCTIONS 

We now turn to the ground state properties of Eq. (JTJ 
as a function of magnetization, concentrating on several 
correlation functions in order to identify the leading in- 
stabilities. Note that our analysis is only valid if M 7^ 0. 
Apart from a term representing the magnetization M in- 
duced by the external field, the longitudinal correlation 
function shows an algebraic decay with distance r: 



M 



d cos(2fc F r + (a - 0)k F ) Kj, 



2n 2 r K + 



87r 2 r 2 
(7) 

The constants Ci, i — 1,2,3, appearing here and in 
Eq. © will be determined through a comparison with 
numerical results. 

In contrast to Eq. ([7]), the transverse zy-correlation 
functions decay exponentially reflecting the gapped na- 
ture of the single spin- flip excitations. Here we do not 



restrict ourselves to the equal-time expression only, be- 
cause we will need non-equal time correlation functions 
to extract the finite-size corrections to the gap later on. 
We obtain: 



<5+(0,0)S^(r,r)) 



_ 1 y e -A 1 (M) y /T^+r^/vl 



(r 2 +v\t 2 Y k + (r 2 +v 2 _T 2 y 



(8) 

where r stands for the Euclidean time, Ai(M) is the 
AS Z = 1 gap, and v± ~ v ± Ji/tt in the weak cou- 
pling limit. The Kronecker delta strictly applies to the 
thermodynamic limit, while on the lattice an additional 
contribution for a 7^ (3 exists. 

It is noteworthy that, different from Eq. ([8]), the in- 
plane correlation functions involving bilinear spin combi- 
nations decay algebraically. This stems from the gapless 
nature of AS Z = 2 excitations. In fact, these are the 
slowest decaying correlators close to the saturation mag- 
netization: 



(5+(r)5+(r)Sr(0)5 2 -(0)) 



C 2 C 3 cos(2fc f r) 

,1/K+ + r K++l/K+ 



(9) 



This result is reminiscent of a partially ordered state be- 
cause the ordering tendencies in this correlation function 
are more pronounced than those of the corresponding 
single-spin correlation function Eq. ([8]). Therefore, we 
call the correlator © 'nematic'. Furthermore, we will 
refer to a situation where Eq. is the slowest decay- 
ing one among all correlation functions as a 'nematic-likc 
phase'. 

By virtue of the exponential decay in ((5|), the correlator 
© is proportional to: 



<(#(r) + #(r)) 2 (Sr(0) + ST(0)) 2 



(10) 



The term (Sf + S%) 2 appearing in the case of the S = | 
zig-zag ladder corresponds to the operator (S a ) 2 in the 
case of a S — 1 chain. One can think of an effective 5=1 
spin formed from two neighboring S = 5 spins coupled 
by the ferromagnetic interaction. A similar behavior of 
correlation functions, namely the exponential decay of in- 
plane spin components and the algebraic decay of their 
bilinear combinations, is encountered also in the AY"2 
phase of the anisotropic 5 = 1 chaini^ and in the spin-1 
chain with biquadratic interactions, see, e.g., Ref. I la. 

The algebraic decay of the nematic correlator as op- 
posed to the exponential decay of (JSJ suggests that there 
are tendencies towards nematic ordering in this phase. 
Depending on the value of K + the dominant instabilities 
are either spin-density-wave ones for K+ < 1 or nematic 
ones for K + > 1 . From the result for K + given in Eq. (J6]) 
one can perturbatively evaluate the crossover value of J\ : 



\Ji, 



1 



irv(m) 
K{m) \K{m) 



1 



(11) 



For J 1 < J\ tCT the nematic correlator © is the slowest 
decaying one, i.e. one is in the nematic-like phase. The 
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FIG. 1: (Color online) Correlation functions at J\ — — J-2 < 0, 
and magnetization M = 3/8: (a) longitudinal component S%, 
(b) transverse component S^r, (c) spin nematic S^rS^r +1 . x is 
the distance in a single-chain notation. ED results for periodic 
boundary conditions are shown by symbols, fits by lines. Note 
the logarithmic scale of the vertical axis in panel (b). 



behavior of the cross-over line can be read off from the 
behavior of K(m): K(m) increases monotonically with 
m, tends to K = 1 for m — > 1/2 , and satisfies K < 1 
for m < 1/2 (see, e.g., Refs. Il6lll7r ). Therefore, we have 
J\,cr — for M = 1/2 with increasing ferromagnetic 
|Ji,cr| f° r decreasing M. This means that for J\ < 
a regime opens at high M where nematic correlations 
given by Eq. ^ dominate over spin-density-wave corre- 
lations given by Eq. ([7]), in agreement with Chubukov's 
prediction^ 

Now we check the correlation functions obtained within 
bosonization against exact diagonalization (ED) results. 
Numerical data obtained for J\ = — J 2 < and M = 3/8 
on finite systems with periodic boundary conditions are 
shown in Fig. [T] This parameter set allows for a clear 
test of the above predictions, but represents the generic 
behavior in the phase of two weakly coupled chains. To 
take into account finite-size effects we use the observation 
that for a conformally invariant theory, any power law 



on a plane becomes a power law in the following variable 
defined on a cylinder of circumference L: 



X 7T 

T 



(12) 



First we fit the nematic correlator given by Eq. f9|, 
which from bosonization is expected to be the leading 
instability at high magnetizations. Using the part with 
x > 5 of the L = 64 data shown in Fig. [TJ;, we find 
1/K+ = 0.904 ±0.011, C 2 = 0.143 ± 0.004, and C 3 = 
— 0.326±0.013. Fig.rTJ; shows that all finite-size results for 
the nematic correlator are nicely described by this fit with 
the dependence on L taken into account by substituting 
Eq. (fT2|) for the power laws. Moreover, from K + > 1, we 
see that the system is indeed in the region dominated by 
nematic correlations for M = 3/8 and J\ — — J2. 

Now we turn to the longitudinal correlation function 
which we fit to the bosonization result Eq. ([7]). Since 
most numerical parameters have been determined by the 
previous fit, only one free parameter is left which we de- 
termine from the numerical results of Fig. [Iji for L = 64 
and x > 14 as C\ — 0.060 ± 0.004. Predictions for other 
system sizes are again obtained by substituting Eq. (fl"2"|) 
for the power laws. The agreement in Fig. [1^, is not as 
good as in Fig. [TJ;. However, it improves at larger dis- 
tances x and system sizes L, indicating that corrections 
omitted in Eq. ([7]) are still relevant on the length scales 
considered here. 

Finally, the xy-correlation function is shown in Fig.[T]3 
with a logarithmic scale of the vertical axis of this panel. 
The exponential decay predicted by Eq. (J5J is verified. 
One further observes that correlations between the in- 
plane spin-operators belonging to different chains (odd 
x) are an order of magnitude smaller than on the same 
chain (even x). This suppression of correlations between 
different chains corresponds to the 5 symbol in ([8]) , which 
strictly applies only in the thermodynamic limit and for 
large distances. 

We summarize the main result of this section: in-plane 
spin correlators are exponentially suppressed for any fi- 
nite value of the magnetization in the parameter region 
\J\\ < Ji. The ground state crosses over from a spin- 
density-wave dominated to a nematic-like phase with in- 
creasing magnetic field, with the crossover line given by 

Eq. cm. 



IV. EXCITATIONS 

We next address the excitation spectrum. Since the 
gap to AS Z — 1 excitations should be directly accessible 
to microscopic experimental probes such as inelastic neu- 
tron scattering or nuclear magnetic resonance, we analyze 
its behavior as a function of magnetization. Sufficiently 
below the fully polarized state the gap can be calculated 
analytically using results from sine-Gordon theory. In 
addition, to leading order of the interchain coupling, one 



can get qualitative expressions using dimensional argu- 
ments for the perturbed conformally invariant model: 



Ai(m) 



c 2 (m)| Ji| sin(7rm) 



v(m)(l — J\K(rn)/'Kv(rn)) 



(13) 



where v = 2 - 2K(m) (l + JiK(m)/Trv(m)) . m(h), K(h) 
and v(h) can be determined numerically from the Bethe 
ansatz integral equationsj 16 ' 17 ! 18 ' 19 

With this information and Eqs. Q and (|13[) we de- 
termine the qualitative behavior of the single-spin gap 
Ai(M) as a function of M: it increases from zero at 
zero magnetization, reaches a maximum at intermediate 
magnetization values, then shows a minimum and, upon 
approaching the saturation magnetization, it increases 
again. As our formulas do not strictly apply at m = 0, 
the notion of a vanishing gap at zero magnetization may 
be a spurious result. Note that when the fully polarized 
state is approached, the magnetization increases in an un- 
physical fashion since in this limit bosonization becomes 
inapplicable. At the point where the magnetization sat- 
urates the exact value of the gap can be obtained from 
the following mapping to hard-core bosons ^ 11 ' 20 



QZ 1 t 



ST = 



(14) 



Comparing Eq. (fl4|) with Eq. (|2|) one recognizes the lead- 
ing terms in Haldane's harmonic fluid transformation for 
bosons. 20 Using a ladder approximation which is exact in 
the two-magnon subspace we arrive at: 



Ai M 



1 



4J 2 



2J1J2 — J 2 J 2 + 8J1 J2 



16J, 2 



2(J 2 - Ji) 
1 Jl(ZJ 2 + Ji) 
8 MJi-Jx) 



8J 2 



(15) 



In Eq. (|15p we have represented the gap as a difference 
of two terms: the quantum and the classical instability 
fields emphasizing its quantum origin. 

In order to verify these field theory predictions, we 
perform complementary numerical computations using 
the DMRG method. 21 Open boundary conditions are im- 
posed and we typically keep up to 400 DMRG states. 
From DMRG we obtain the ground-state energies E(S Z ) 
as a function of total S z . For those values of S z that 
emerge as a ground state in an external magnetic field 
we compute the single-spin excitation gap from 



Ai(M) 



E{S Z + 1) + E(S Z - 1) - 2 E(S Z ) 



(16) 



Fig. [2^ shows numerical results for Ai at a selected value 
of Ji = —0.3 J2 < for the largest system sizes investi- 
gated. We find that the finite-size behavior of the gap 
Ai(M, L) for system sizes L > 24 is well described by 
a 1/L correction. This will be further corroborated by 
field-theoretical arguments outlined below. Therefore, we 
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FIG. 2: (Color online) Density matrix renormalization group 
results for the gaps at Ji = —0.3 J2 < as a function of 
magnetization M. Panel (a) shows the single-spin excitation 
gap ()16p . panel (b) the finite-size gap (|19p for two flipped 
spins multiplied by the chain length L. 



extrapolate it to the thermodynamic limit using a fit to 
the form 



Ai(M,i) = Ai(M) 



a(M) 



(17) 



allowing for an additional 1/L 2 correction for those values 
of M where at least 4 different system sizes are available. 

This extrapolation is represented by the full circles in 
Fig. [2^; errors are estimated not to exceed the size of the 
symbols. Our extrapolation for Ax is consistent with a 
vanishing gap at M — in agreement with previous nu- 
merical studies^ 2 - although bosonization predicts a non- 
zero - possibly very small - gap i 13 i 22 ' 23 The behavior of 
Ai(M) confirms the picture described above: the gap is 
non-zero for M > 0, goes first through a maximum and 
then a minimum and finally approaches A1/J2 ~ 0.023 
given by Eq. for M -> 1/2. 

We further wish to point out that for chains with pe- 
riodic boundary conditions, the coefficient a(M) of the 
finite-size extrapolation Eq. (| 17[) is determined by the 
spin-wave velocity and the critical exponent of the soft 
mode from the AS Z = 2 channel. Indeed, using Eq. |8]) 
where we can set r — 0, and use the conformal mapping 
(fT2")l to the cylinder, we see that the leading finite-size 
correction to the gap is: 



Ai(M,Z)=Ai(M) + - 



1 irv + (M) 



L4K + (M) 



(18) 
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FIG. 3: (Color online) Numerical dispersion spectrum in the 
subspaces of odd S z computed for L = 24 and Ji = — J2 < 0. 
The wave vector is given relative to the ground state wave 
vector (0 for S z = 0, 4, 8 and tt for S z = 2, 6). 



Note that we have to replace sin with sinh in Eq. (Q2 
in order to extract a gap, since we are dealing with Eu- 
clidean time. In addition we used the fact that in our ap- 
proximation the effective Hamiltonian ([5]) is a direct sum 
of symmetric and antisymmetric sectors. Moreover, it is 
only the symmetric sector enjoying conformal invariance 
and consequently we perform the replacement t — * sinh r 
only in the symmetric sector. The antisymmetric sector 
has a spectral gap and its contribution to the finite-size 
corrections of the single-spin flip excitation energy are 
exponentially suppressed with system size<21 With this 
method one cannot fix the amplitudes of the 1/L 2 term 
and beyond. Note furthermore that there may be ad- 
ditional surface terms for open boundary conditions as 
employed in the numerical DMRG computations. Nev- 
ertheless there is a dominant 1/L correction in any case. 

Next, we briefly look at the AS Z — 2 excitations. Their 
finite-size gap is, in analogy to Eq. (|16|) . computed with 
DMRG from 



A 2 (Af) 



E(S Z 



E{S Z 



2E{S Z 



(19) 



Fig. [2}d shows numerical results for L A%{M, L) again at 
the value Ji = —0.3 J 2 < 0. One observes that the scaled 
finite-size gaps collapse onto a single curve which shows 
that A2(M, L) scales linearly to zero with 1/L, exactly 
as expected for gapless excitations in ID. Furthermore, 
we observe that the scaled quantity L A%(M, L) vanishes 
as one approaches saturation M = 1/2 which indicates a 
vanishing of the velocity of the corresponding excitations 
at saturation. 

We proceed by discussing the wave-vector dependence 
of the AS Z = 1 excitation, while we remind the reader 
that the low-energy excitations are in the AS Z = 2 sec- 



tor. Fig. [3] shows representative ED results obtained for 
rings with L = 24 and J\ = —J2 < 0. For ground 
states with low S z , the AS Z = 1 excitation spectrum 
looks similar to the continuum of spinons. On the other 
hand, close to saturation one has single-magnon excita- 
tions with a minimum given by the classical value of 
the wave vector k c i = arccos(| J\ |/4 J 2 )£ i 13 i 25 We read 
off from Fig. [3] that upon lowering the magnetic field, 
this minimum shifts from the classical incommensurate 
value towards 7r/2, i.e. the value appropriate for two de- 
coupled chains. This renormalization of the minimum of 
the magnon excitations towards the value of decoupled 
chains can be interpreted in terms of quantum fluctua- 
tions, which are enhanced when the density of magnons 
increases. A strong quantum renormalization of the pitch 
angle from its classical value at zero magnetization was 
previously observed by the coupled-cluster method and 
DMRG calculations^ 



V. SUMMARY 

We have combined numerical techniques with analyti- 
cal approaches and mapped out the ground state phase 
diagram of the frustrated ferromagnetic spin chain in an 
external magnetic field. We have established that with 
increasing magnetic field, the ground state crosses over 
from a spin-density-wave dominated to a nematic-like 
phase. Single spin flip excitations are gapped, giving rise 
to an exponential decay of in-plane spin correlation func- 
tions in both regimes. We have studied the single- and 
two-spin flip excitation energy numerically. Using tools 
from conformal field theory we have further shown that 
the amplitude of the leading 1 / L correction term to the 
single-spin Hip gap is determined by the critical exponent 
and the spin-wave velocity of the soft mode. 

Finally, in order to apply our findings to the mate- 
rial LiCuV04, one should take into account interchain 
interactions as well as anisotropies, which are expected 
to be present in this system.— At low fields, a helical state 
has been observed experimentally.-^ On the other hand, 
for the purely one dimensional case, we have shown that 
upon increasing the magnetic field there is a competition 
between spin-density-wave and nematic-like tendencies. 
Those are the two leading instabilities at high magnetiza- 
tions and thus they are the natural candidates to become 
long-range ordered in higher dimensions. The question 
whether there are true phase transitions at high fields 
in higher dimensions is beyond the scope of the current 
work. 
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